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How statistical forces depend on thermodynamics and kinetics of driven media 

Urna Basu, Christian Maes 
Instituut voor Theoretische Fysica, KU Leuven, Belgium 

Karel NetocnjQ 

Institute of Physics, Academy of Sciences of the Czech Republic, Prague, Czech Republic 

We study the statistical force of a nonequilibrium environment on a quasi-static probe. In the 
linear regime the isothermal work on the probe equals the excess work for the medium to relax to 
its new steady condition with displaced probe. Also the relative importance of reaction paths can 
be measured via statistical forces, and from second order onwards the force on the probe reveals 
information about nonequilibrium changes in the reactivity of the medium. We also show that 
statistical forces for nonequilibrium media are generally nonadditive, in contrast with the equilibrium 
situation. Both the presence of non-thermodynamic corrections to the forces and their nonadditivity 
put serious constraints on any formulation of nonequilibrium steady state thermodynamics. 


Statistical forces appear in macroscopic physics as gra¬ 
dients of thermodynamic potentials [I]. The equilibrium 
pressure of a gas in a container is for example given in 
terms of the change in free energy under a small exten¬ 
sion of the volume, an essential tool for obtaining the 
equation of state. In irreversible thermodynamics such 
statistical forces translate the typical tendency of macro¬ 
scopic systems to reach equilibrium and generate fluxes 
of matter, momentum or energy [2j. Another simple and 
speaking example is that of the entropic force in a spring 
as governs basic polymer physics [3j. 

It is of much current interest to extend the study of 
statistical forces to nonequilibrium regimes EH]. There 
are multiple motivations beyond the general ambition of 
constructing a nonequilibrium statistical mechanics. Un¬ 
derstanding the nature of pressure in active media or the 
change in osmotic pressure in presence of active solutes is 
fundamental for the viability of steady state thermody¬ 
namics HH9j. There are also a growing number of practi¬ 
cal aspects of using nonequilibrium environments to mod¬ 
ify the nature of statistical forces da- For example, it 
is important to explore the possible phenomenology of 
Casimir forces m or of forces derived from vector and 
higher order potentials that induce oscillatory behavior 
on the probe’s motion. Finally, and as we show, statisti¬ 
cal forces give better operational meaning to the notion 
of dynamical activity as recently explored for nonequilib- 
rium purposes mm- 

This Letter is a systematic theoretical account on the 
meaning of statistical forces within the standard set-up of 
classical mesoscopic systems. Our main findings are first 
that the linear order correction to the equilibrium statis¬ 
tical force is still of thermodynamic nature (but not in 
the form of a gradient). In that regime statistical forces 
measure excess work which is a key-quantity in present 
accounts of steady-state thermodynamics 0 m ■ Sec¬ 
ondly, we find that statistical forces can measure dynam¬ 
ical activity. That includes both measuring the relative 
activity through different reaction paths, and specifying 
how the external driving affects reactivities. Thirdly we 


discuss why statistical forces outside equilibrium become 
nonadditive. That connects with typical long range ef¬ 
fects in nonequilibrium. 

All of these aspects yield basic understanding of statis¬ 
tical forces outside equilibrium and can already be illus¬ 
trated with toy models of mesoscopic systems that catch 
the essential physics. 


FIG. 1. (Color online) (a) A slow probe (light gray disc) is 
immersed in a nonequilibrium medium (dark green circles), 
in contact with an equilibrium reservoir (small blue circles), 
(b) A dilute colloid solution diffusing on a ring, each colloid 
particle being coupled to an internal degree of freedom which 
is driven out of equilibrium. 

Set-up: Imagine a nonequilibrium medium with a family 
of driven degrees of freedom. It for example consists 
of particles undergoing nonconservative forces dissipating 
into a thermal bath at inverse temperature /3. A slow 
probe with position x (a possibly high-dimensional or 
even collective coordinate) is immersed in the medium, 
with interaction potential U{x,rf), which is assumed also 
to include the self-interaction of the medium. Fig. [l|a) 
shows a representative cartoon of the set-up. We study 
the statistical force 

f{x) = - J Pxid'q) V x U{x,r]) = -iy x U(x,rj)) x 

where the average is over the steady nonequilibrium sta¬ 
tistical distribution p x of the ^—medium at fixed x. We 
assume there that the variables ?y are relaxing very fast on 
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the typical time-scale of motion of the probe position. In 
other words, the probe moves in a medium which is in in¬ 
stantaneous stationarity, be it nonequilibrium. Note that 
the above formalism employs Newton’s action-reaction 
law between probe and medium. 

When the medium dynamics satisfies detailed balance, 
the statistical force / eq ( x) = ~'V x T(x) is given in terms 
of the free energy T{x) = — log Z x with the equilibrium 
partition function Z x corresponding to the 77-medium 
when in equilibrium with fixed probe position x. The 
ensuing motion of an overdamped probe would be gra¬ 
dient flow along x = — S 7 x F{x) (on the appropriate slow 
time scale) to reach equilibrium at the position x that 
minimizes the free energy. 

Before discussing the thermodynamic and kinetic 
aspects of statistical force from a driven medium in 
more detail we illustrate the set-up and issues with 
the help of a simple example of a diffusion coupled to 
nonequilibrium degrees of freedom. 

Example — diffusion coupled to rotator —The probe is a 
colloid trapped on the unit circle, x £ S 1 , in contact with 
a thermal bath at inverse temperature (3 = 1 , modeled 
via the overdamped Langevin dynamics 

7i:= + ]f^ ^ W 

where 7 is the damping and is standard white noise. It 

is coupled to an ‘internal’ degree of freedom 77 = —1,0,1 
(Fig. [ljb)) through the potential 

U(x,ri) = ijsinx+ 2 p 2 cos x ( 2 ) 

The choice of potential is rather arbitrary but avoids spe¬ 
cial symmetries. The fast degree of freedom 77 is driven 
and rotates with transition rates 

k x {r],r]') = ( 3 ) 

The drive affects both the antisymmetric 5(77,77') = 
—s(rf,T)) and symmetric 0(77,77') = 0(77', 77) parts of 
the rates. We consider a uniform drive with field e, 
s(—1,1) = s(l,0) = s(0, —1) = / 3 s and for simplicity 
we assume that only one of the reactivities gets changed 
by the drive 0(—1,1) = 0(1, —1) = 0o(l + a|e|) for some 
a > 0 , while 0 ( 0 , ± 1 ) = 0 (± 1 ,O) = 1 . The equilibrium 
(detailed balance) reference has e = 0 . There, the 0 O 
picks up the relative importance of the dynamical activ¬ 
ity over the transitions lf>—1. 

The statistical force on the particle is 

f(x) = — (77)* cosx + 2 {r] 2 ) x sinx ( 4 ) 

The equilibrium force f e q{x ) and nonequilibrium correc¬ 
tion g(x) = f(x) — f e q (x) are plotted in Fig. [ 2 ](a). While 
the force is derived from a potential in equilibrium that 
is no longer true when the 77 are driven; the rotational 
part of the force / ro t = / f(x) da; is plotted versus e in 
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FIG. 2. (Color online) (a) The equilibrium statistical force 
feq{x) and the nonequilibrium correction to it g(x) for field 
s = 2.0. Here 0o = l,a = 2. The inset shows the rotational 
part of the force f IO t as a function of the field e. (b) The 
excess work dW e *(x) is the extra dissipated work for relaxing 
to the new steady nonequilibrium (at different power W) as 
the probe position x —> x + dx has changed. 

the inset. 

Perturbative approach: We continue with the general 
set-up. The equilibrium distribution p eq (drf) satisfies the 
Gibbs formula 

V^ q (d77) = —p eq (drj) [f 3 S 7 x U(x,rj) +V*log Z x ] ( 5 ) 

Multiplying the above relation with p x (dr])/p% q (dr)) : 
where p x (drj) is the stationary nonequilibrium den¬ 
sity and integrating over 77 we find the force f(x) = 
j]X 7 x log Z x + g(x) with nonequilibrium correction 

«(.)—i/«■(*,) (6) 

as an expectation with respect to the equilibrium refer¬ 
ence. The additional nonequilibrium force hence derives 
from the nonequilibrium correction h x in the stationary 
density, p x (drj) = p% q (drj) [1 + h x {rjj\, as does the work 
(scalar product of force with displacement) performed on 
the probe, 

g(x) ■ dx = J p eq { d??) h x+dx {rj) ( 7 ) 

No close-to-equilibrium assumptions have been made so 
far. Yet, © 0 invite building perturbation expansions 
g{x) = eg\{x) +s 2 g2(x) + ... around the equilibrium case 
s = 0 . For such an expansion we assume that the medium 
satisfies the local detailed balance condition d El, 
which allows to quantify the amount of nonequilibrium 
driving in terms of a parameter s and through which the 
(extra) entropy flux to the thermal bath can physically 
be identified as for example the 5(77,77') in ([ 3 |. 

The first-order nonequilibrium correction in the sta¬ 
tionary distribution is given by the McLennan for¬ 
mula dug, which allows a thermodynamic character¬ 
ization of h x {rf). More precisely, always in the linear 
regime, h x /f 3 is given in terms of the excess work as has 
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appeared in constructions of steady state thermodynam¬ 
ics mm- The notion of excess work arises from the fact 
that the nonequilibrium medium for fixed x constantly 
dissipates work W (as Joule heating) even at stationar- 
ity. When now the probe position x —> x + da; changes, 
the medium starts to relax to a new stationary condition 
Px —► Px+dx with extra, indeed excess, work 3 W e *{x) 
done by the non-conservative forces; see Fig §b) for a 
representation. That is of course also dissipated as heat 
in the thermal environment. The McLennan scheme then 
applies to ([ 7 ]) and readily implies that the work done on 
the probe by the medium equals the excess work done on 
the medium by the nonequilibrium drive, 

9 l (x) ■ dx = dW ex {x) (8) 

That identity is the first result of the present Letter; 
it is not true in general beyond the linear regime. The 
result can be extended in the same regime of weak 
nonequilibrium thermodynamics to include processes 
that are not isothermal using quasi-static energetics. 


Kinetic aspects: As a new application of statistical 
forces we find that they also give information about the 
dynamical activity in the medium. Just like excess work 
above, notions of dynamical activity while clearly rele¬ 
vant for nonequilibrium studies so far lack operational 
meaning. 

The first point is that we can infer from the linear order 
g i information on the equilibrium reactivities. We show 
it with the example in ([ 3 ]). The relative equilibrium fre¬ 
quency of undirected transitions between states is there 
governed via the <f> o- While that <f>o is not influencing the 
equilibrium free energy, it does become visible to first or¬ 
der gi in the statistical force. In other words, the relative 
dynamical activity in paths of equilibrium media can be 
found from applying a small force and from measuring 
the resulting mechanical effect on the probe. Fig. [ 3 ]) a) 
shows the dependence on 0 0 in the slope of the statisti¬ 
cal force g(x = 7t/2) at small e; r(cf> o) = gi(n/ 2 ). The 
more general dependence is 


9i{x) 


A x + B x (f> o 
C x + (f> o 


where the A,B,C are only depending on x (and not on 
a for example). 

But there is more: the second order is able to pick 
up the e—dependence (parameter a) in the reactivities, 
invisible to linear order. That is in line with the analysis 
of higher order effects in the response formalism in d 3 
The linear and second order corrections in the statistical 
force generated by ([ 3 ]) are plotted separately in Fig. [ 3 ])b). 
The statistical force starts to depend on the change a 
in reactivity only from second order onwards around 
equilibrium. Alternatively, from the force on the walker 
we can detect information about the reactivity-change 



FIG. 3. (Color online) (a) The statistical force g(x = 7t/2) 
as function of e. The initial slope depends on the relative 
equilibrium reactivity 0o (here 0 o = 0.5). For larger drive e 
the dependence on the reactivity parameter a becomes vis¬ 
ible. (b) Linear gi(x) and second order correction g 2 {x) as 
a function of the probe position for a = 1,2 and 0 o = 1 . 
The linear order is independent of a while the second order 
depends explicitly on it. 


under the driving field. 

Nonadditivity: The medium rj can consist of multi¬ 
ple reservoirs which are spatially separated or are only 
in weak contact. We know that in the equilibrium case, 
because of the locality of the interaction or from the ad¬ 
ditivity of the free energy T{x), the statistical forces add 
the total force on the probe is the sum of the forces 
from the different reservoirs. To be specific let us think 
of the situation where the probe is coupled to two dif¬ 
ferent media both of which can be driven out of equi¬ 
librium in general. From equation © , we see what is 
needed for additivity of statistical forces: the station¬ 
ary 77—distribution must approximately factorize over its 
constituents. Or, the coupling with the second reservoir 
should not completely change the stationary distribution 
in the first reservoir. That is exactly what does not need 
to be the case for coupled nonequilibrium reservoirs, and 
thus can lead to nonadditivity of forces. 



FIG. 4. (Color online) (a) Coupling of two internal degrees 
of freedom to the probe. The A-reservoir is driven and will 
create a current also in the B-reservoir, making the forces 
nonadditive, (b) Nonadditivity in the rotational component 
of the total force A/ rot versus e for 7 = 1.0, 5.0 and a = 
2,</>o = 1. In the inset: the deviation from additivity Af(x) 
versus the probe position x at e = 2 . 0 . 

In the simplest case we have the probe in contact with 
independent particles distributed over two parts, sub- 
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medium A and submedium B. The total statistical force 
is then of the form 


f = PA Ia + Pb Sb ( 9 ) 

where pa,b are the concentrations and Ja,b are the forces 
per particle. There is a parameter 7 > 0 that couples 
A and B, and there is a nonequilibrium driving e that 
either drives the particles in A, or in B, or in both parts. 
Particles can move between A and B as long as 7 > 0 . 
Both the concentrations p e AB and the forces per particle 
/as generally depend on the driving e as well as on the 
coupling 7. If the connection gets interrupted, setting 
suddenly 7 = 0, then we generally observe a shift in 
the statistical force, even if the concentrations remain 
unchanged. Indeed, by interrupting the flow between A 
and B the currents and the stationary distributions may 
change in both submedia, and the forces per particle take 
on new stationary values. The observed change in the 
force or measure of nonadditivity is then defined as 

A / = pT [.fT - f e A°] + Pb 1 [/s’ 7 - /s’ 0 ] ( 10 ) 

To illustrate the nonadditivity of nonequilibrium sta¬ 
tistical forces we take again the example of a diffusive 
probe, now coupled to two reservoirs rj = {t)a, Pb}- 
For simplicity we again consider each of the reservoirs 
as a three state rotator with 77 = —1,0,1 as before. 
The probe is coupled to the media through energy func¬ 
tion U(x,rj a ) = 5 a ^ A [r]SV[vx + 2 rj 2 cosx] + 8 a ^ B [ricosx + 
2 rj 2 sin x] which is a sum of potentials similar to ([2]) . 
Additionally the two reservoirs are also coupled to each 
other via two distinct ‘bridges’ at 1a -H- Is and — 1a <-)• 
— 1 b with jump rates 

k x (±l A , ± 1 B ) = 7 [u{x,±i B )-u{x.±i A )} 

and similarly for B —> A; see Fig. B a ) which repre¬ 
sents a doubling of internal degrees of freedom with re¬ 
spect to the situation depicted in Fig. [IJb). The whole 
system can be thought of as the probe being coupled 
to indistinguishable independent walkers hopping on the 
6 states p a with 77 = — 1 , 0,1 and a = A,B. Writ¬ 
ing Za,b = SijeA s the equilibrium statistical 

force of the form | 9 | equals 

fe q = - p V, log(^ + Z B ) = p° A f° A + p° B f°B 

where f AB = jj log Z A B is the mean force per particle 
from A and B, respectively, and p\\ B = Z A B /(Z A + Z B ) 
is the probability of the particle to be in A or B. In equi¬ 
librium the statistical force is strictly additive because 
under detailed balance the prefactor 7 does not change 
the Boltzmann occupation statistics; the f A B do not de¬ 
pend on the coupling 7 so that by cutting the connection 
between A and B nothing is changed. However, when 
either A or B (or both) are driven out of equilibrium 


the bridges and the transition strength 7 become impor¬ 
tant. For example, a drive in medium A typically induces 
a current in B and changes the stationary distribution, 
thus giving rise to a different statistical force from B. The 
nonadditivity ( 10 ) is shown for that example in Fig. [ 4 ](b) . 
The rotational component of the difference A/ is plotted 
there for different 7. The difference A f(x) for a fixed 
driving field e in A as in ([ 3 | is shown in the inset. As 
expected | A/| increases with the strength of the coupling 
7 - 


The fact that nonadditivity already occurs for nonin¬ 
teracting particles is perhaps surprising. It is however 
directly related to the phenomenon that local changes in 
systems with conserved quantities can have global effects 
on the nonequilibrium stationary distribution, [21] . 

For interacting particles in the media with generic 
long range correlations, that effect can only be enhanced 
and other issues become also important. Indeed, as 
a related effect for spatially extended systems, the 
difference between bulk and boundary contributions to 
statistical forces is fading [ 3 ]. All that is relevant for the 
status and possibility of steady state thermodynamics, 
the macroscopic theory of nonequilibrium systems. In 
particular, the difference between bulk versus surface 
effects is important in the usual thermodynamics, and 
the equations of state typically connect pressure with 
bulk properties. The above simple arguments and exam¬ 
ple already indicate possible limitations for macroscopic 
nonequilibria. 


Conclusion: Statistical forces on quasi-static probes in 
nonequilibrium media have a number of new interesting 
features with respect to the equilibrium case. In lin¬ 
ear order around equilibrium they can still be obtained 
from irreversible thermodynamics, where the force mea¬ 
sures the excess work of driving forces on the medium. 
Note that already there the forces can be rotational and 
nonadditive as a result of long range effects. For non¬ 
thermodynamic information these forces can be used to 
measure relative and excess dynamical activities in the 
medium, a concept of growing importance for which lit¬ 
tle operational tools are available so far. Nonequilibrium 
conditions might affect reactivities and time-symmetric 
aspects of the medium, and they are present in the statis¬ 
tical force on the probe from second order onwards. Such 
higher order effects are due to the coupling of the en¬ 
tropy flux with excess dynamical activity, [ 12 | . That also 
contributes to the violation of the second fluctuation- 
dissipation relation in generalized Langevin equations for 
a probe coupled to nonequilibrium media, [22]: there 
again the friction picks up more kinetic information from 
the medium. Statistical forces thus provide an important 
tool to probe the nature and properties of nonequilibrium 
systems. 
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